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Abstract 

The angular and the radial parts of the dynamics of the perturbed Kepler motion are 
separable in many important cases. In this paper we study the radial motion and its 
parametrizations. We develop in detail a generalized eccentric anomaly parametriza- 
tion and a procedure of computing a generic class of integrals based on the residue 
theorem. We apply the formalism to determine various contributions to the luminosity 
of a compact binary. 



1 Introduction 



The main sources of the gravitational radiation to be detected by Earth-based 
observatories (LIGO, VIRGO, TAMA and GEO) are compact binaries. The 
Laser Interferometer Space Antenna (LISA) will also detect gravitational waves 
from colliding galactic black holes. Compact binaries are composed of black 
holes and/or neutron stars. In such syst ems, a h ighly ac curate descri ption would 
include the spin-spin (SS) interaction [Barker & O'Connell I ll 1 979l)l . the mag- 
netic dipole-magnetic dipole moments (DP) co ntribution lloka fc Taniguchi I 
()2000|) ] and the quadrupole-monopole effect (QM) IPoisson I jl998|)l . The parametriza- 
tion o f the ra dial motion for them are prese nted in 
fcOOSl) ]. and jOergelv fc Keresztes I fefUlfl ]. 
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The method for sol ving: a wide class of the perturb ed Keplerian radial mo- 
tions is worked out in Gergelv. Peries. Vasuth I (j2000)]. For this purpose two 
(the true and eccentric anomaly) parametrization were introduced. The method 
can be further generalized for an even wider class of perturbed Keplerian radial 
motions (see: Gergelv. Keresztes. Mikoczi 1 1|2006|) ]'). 



2 The perturbed Kepler motion 

The so-called radial equations for S S, DD and QM interactions are given in 
[Keresztes. Mikoczi . Gergelv I l)2005|) ]. T he radial equation of the generalize d 
perturbed Kepler motion we consider here [Gergelv. Keresztes. Mikoczil l|2f)f)fil) ] 
is 

2E , 2Gm L 2 , A <^(x) 

2^ ,,2 r i y 1 ) 



r 2 



fi r /i 2 r 2 [i^r 



where fi(x) characterize the perturbative terms and they are periodic functions 
of the true anomaly: 



<P. 



(x) = (A/ + 9ij cos x) sin J X ■ (2) 



3=0 



The expression (Pi(x) given a bove is equivalent with a generic Fourier series 
[Gergelv. Keresztes. Mikoczil l[200fi|) ]. The coefficients and #y can be ex- 
pressed in terms of the coefficients of the Fourier expansion as we ll. The last 
term in Eq. Q contains the generic perturbing Br umberg force iBrumberg 



lll99lln. the spin-orbit interactio n for compact binaries [Rieth and Schafer Il[l997 ) 



iGergelv. Peries. Vasuthl l[l998h ]. and the SS, DD and QM contributions. The 
energy E and angular momentum L refer to the perturbed motion. From the 
condition r 2 — we find the turning points r max : 



Gmfj,±A 1 ±/ \ 

-2E ± ^A ^ {X) 



L 2 



(3) 



1 /2 

where Aq — (G 2 m 2 /i 2 + 2EL 2 //j) is the magnitude of the Laplace-Runge- 
Lenz vector belonging to the perturbed motion characterized by E and L. ip^ = 
<Pi(0), (fit = tpiift) are small coefficients. With these turning points is possible 
to introduce the generalized true anomaly parametrization of the radial motion, 
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• (x) as 

2 1 + cos x 1 - cos x , . s 

- = H . (4) 

' ' m i n ax 



Integrals of the type 



„2+n 



dt (5) 



frequently occur. Here w(x) is the same type of periodic function of the true 
anomaly as (fi (%), and T is the radial period of the motion. 

For constant values of ifi and n > these integrals can be evaluated in terms 
of a complex true anomaly variable z = exp (iy) by using the residue theorem. 
The only pole is in the origin (z = 0) Gergelv. Peries. Vasuth I l|2fl0flh ]. For the 



class n < a complex eccentric anomaly parameter can be used and then the 
poles are in the origin and in 

1/2 



[ Gmn 2 - yZ-2^ 2 \ 

w\ = —, - , (6) 

\Grnu, 2 + ^/-2u.EL 2 ) 

howev er the latter occurs only rarely in physical applications ( see: [Gergelv. Peries. Vasuth I 

feoorti p. 



3 Generalized perturbed Kepler motion with the 
eccentric anomaly parametrization 

We introduce the r (£) eccentric a nomaly parametrization in the same way as 
in [Gergelv. Peries. Vasiithl l(2f)0(Tl) ]: 

2r = (1 + cos£)r min + (1 - cos£)r max . (7) 

We use the relations between the true and eccentric anomaly |@J, Q 



Gmfj, cos £-A . V^f^ sin ^ 

cos ^ = 7=; a 1' Sm x = -Fr - A 7 ■ ( 8 ) 

um/j — Aq cos 4 Gm/i — Ao cos £ 

Thus we can express (fi as the function of £: 

j=£ V Gm/i - cos J V Gmp? - A u, cos J ' [ ' 
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Employing the eccentric anomaly parametrization l|7jl. the integrals (0 could 
be evaluated as 

r27r (Idt 



„n+l 



For nl = — n — 1 > we apply the binomial expansion 



(10) 



(11) 



leading to a polynomial in cos £. From the radial equation using the eccentric 
anomaly parametrization JZJ to leading order we obtain: 



1 dt 

rd£ 



-2E 



E 



2Mo sin 2 £^ 



(?) 



ri-2 



(12) 



with 



tt s ± = (-^) [^+(Gm M -Ao) t - 2 ±^(Gm Ai + J 4o) i 



(13) 



The bracket becomes proportional to sin 2 £ in (|12(l if the following two conditions 
are satisfied: 

J2 J ^(Gm f x±A ) i (f a ±ga)=0. 



i=0 



(14) 



With conditions l|14|) satisfied, using the parametrization (JJJ, the integrand 
of (|10|l becomes regular. The conditions are fulfilled in the case of the SS, DD 
and QM interactions. Introducing the complex variable w = exp (z£), the inte- 
gral of (|10|) contains only two poles: the origin and the w\ (see: Gergelv. Peries. Vasuth I 
( 2000)]). We have proven for the n < case: 

Theorem: For all perturbed Kepler motions characterized by the radial 
equation with periodic perturbing functions </?j(£) obeying the conditions 
i|14|). and for arbitrary periodic functio ns q>(£), the integrals \ l( Jb are g iven by 
the residue in the origin and in the w\ \Gereelv. Peries. Vasuth I \ 



|)J. on 
complex plane. 

In the next section we apply the above method for computing different con- 
tributions to the luminosity of compact binaries. 
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4 Application for some compact binaries 

Peters and Mathews jPeters and Mathews I l|l963h ] computed the luminosity 
of the compact binaries in the Kepler motion (Newtonian case) with orbital 
parameters: 

G 4 mV 



Cjst — 



15c 5 a 5 (l - e 2 ) 7 /'- 



(37e 4 + 292e 2 



96 



(15) 



where a the semi-major axis and e the eccentricity. The form of ifi (£) can be 
derived for all of the SS, QM and DD contributions. We can then compute the 
various contributions to the luminosity (— {dE/dt}, the time-averaged energy 
loss over one radial period). These are: 

G A m 2 nS\Si 



C 



SiSi 



480c 7 a 7 (l - e 2 ) 10 / 2 



c\ sin K\ sin k 2 cos 2(ipQ — ip) 



C 



QM 



+C2 COS K\ COS K2 + C3 COS 7 

G 4 mV 
60c 5 a 7 (l - e 2 ) 10 /2 



/I Pi [(04(3 sin 2 Ki — 2) + C5 sin 2 Ki cos^j)] 



C 



DD — 



G 3 m 2 ndid 2 
30c 5 a 7 (l - e 2 ) 10 /: 



:[ciAo - c 5 B ] , 



(16) 



where we have denoted by Si the spin magnitudes, by di the magnitudes of the 
magetic dipoles, by pi the quadrupole mom ent scalars and b y .An , Bp, (5,;, 7, ^0, 



ii auxilia ry angular quantities d e fined in [Gergelvl l|200flh ]. |Vasuth. Keresztes. Mihalv. Gergelvl 
])] and [Gergelv fc Keresztes I <|2003h ] . The constants C1..5 are 



2;; 



(17) 



with coefficients Cy given in Tab. ^ 



5 Summary 

We have introduced a generalized eccentric anomaly parametrization for the 
perturbed Kepler motion. We have proved that even for the generic perturbation 
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Table 1: The coefficients in the dj. 








1 


2 


3 


1 





131344 


127888 


7593 


2 


-124864 


-450656 


-215544 


-8532 


3 


42048 


154272 


75528 


3084 


4 





8208 


7988 


474 


5 


2600 


9376 


4479 


-177 



functi ons ipi (£) there are no new poles as compared to Gergelv. Peries. Vasuth I 
(2000)]. The method of integration can be widely employed in the case of 
compact binaries. We have applied the procedure to compute the SS, QM and 
DD contributions to the luminosity of a compact binary. 
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